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Abstract 

The nature of the internal thermal X-ray emission seen in "thermal composite" 
supernova remnants is still uncertain. Chandra observation of the 3C391 shows a 
southeast-northwest elongated morphology and unveils a highly clumpy structure of 
the remnant. Detailed spatially resolved spectral analysis for the small-scale features 
reveals normal metal abundance and uniform temperature for the interior gas. The 
properties of the hot gas comparatively favor the cloudlet evaporation model as a 
main mechanism for the "thermal composite" X-ray appearance, though radiative 
rim and thermal conduction may also be effective. A faint protrusion is found in Si 
and S lines out of the southwest radio border. 

Key words: radiation mechanisms: thermal, supernova remnants: individual: 
3C 391 (G31.9+0.0), X-rays: ISM 
PACS: 95.30.Gv, 98.38.Mz, 95.85.Nv 



1 Introduction 



3C 391 (G31. 9+0.0), with irregular morphology (e.g. Chen & Slane 2001), has 
been classified into the "thermal composite" category of supernova remnants 
(SNRs) (Rho & Petre 1998). They generate bright thermal X-ray emission in- 
terior to their radio shells, and have faint X-ray rims. They are usually found 
to interact with adjacent molecular clouds, characterized by the hydroxyl rad- 
ical maser emission (e.g. Yusef-Zadeh et al. 2003). The nature of the internal 
thermal X-ray emission seen in mixed morphology remnants is still uncertain. 
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So far at least four candidate scenarios compete to account for centrally- 
brightened X-ray morphology. The first scenario is radiative cooling of the 
rim gas. Under this hypothesis, the gas at the rim has been cooled down in 
the radiative stage, with a temperature so low that its X-ray emission is very 
weak, while the gas in the inner volume is still hot enough to emit strong X-rays 
(e.g. Rho & Petre 1998). The second mechanism invokes thermal conduction. 
It is suggested that thermal conduction in the remnant can prevent formation 
of the very tenuous, hot gas in the inner part and therefore change the interior 
structure from the standard Sedov solution, resulting in a nonnegligible density 
and luminous X-ray in the interior (as observed in, especially, the radiative 
stage) (Cox et al. 1999). The third scenario is cloudlet evaporation in the SNR 
interior. When an SNR expands in an inhomogeneous interstellar medium 
(ISM) whose mass is mostly contained in small clouds, the clouds engulfed 
by the blast wave can be evaporated to gradually increase the density of 
the interior gas; consequently, the SNR appears internally X-ray brightened 
(White & Long 1991). The fourth suggestion is that the mixed morphology is 
a projection effect. For shell-like SNRs that evolve in a density gradient such 
as at the edge of a molecular cloud, if the line of sight is essentially aligned 
with the density gradient as well as the magnetic field, the SNRs will appear as 
thermal composites (Petruk 2001). Here we report on high resolution Chandra 
observation of 3C 391 and compare the properties with these scenarios. 




Fig. 1 Smoothed diffuse emission from SNR 3C 391 in the broad band 0.3-7.0 keV (with count-to-noise ratio 
of 6). The color is logarithmically scaled in the range (1.19-331.85) xlO -2 ct s" 1 arcmin . The location of 
the sources removed from the data before the smoothing are marked by white cross labels. All the regions 
used for spectral analysis are indicated in white, with cyan numerical labels. The overlaid 1.5 GHz radio 
contours (in green) are at 1.5, 4.5, 13.6, 28.8, and 50xl0~ 3 Jy beam" 1 (Moffett & Reynolds 1994). The 
two yellow plus signs denote the OH maser points (Frail et al. 1996). The white dashed lines denote the 
border of the S3 chip. 
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2 Data and analysis 




a (2000) 



Fig. 2 Smoothed narrow band 1.2-1.5, 1.7-2.0, and 2.3-2.6 keV (including Mg Hea, Si Hea, and S Hea, 
respectively) diffuse emission images (with S/N ratio of 3) overlaid with the dashed contours of 1.5 GHz 
radio emission (at 1.5, 3.2, 8.25, 16.7, 43.7, and 62.2xl(T 3 Jy beam" 1 ) (Moffctt & Reynolds 1994). The 
seven levels of solid contours arc plotted with square-root intensity scales between the maximum and the 
15% maximum brightness. The two plus signs in each panel denote the OH maser points. 
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Fig. 3 EW images of Mg, Si, and S lines. The Mg and Si (S) images are extracted with 4" (8") pixels and 
smoothed by a Gaussian with a = 12" (16"). The seven levels of solid contours are plotted with square-root 
intensity scales between the maximum and the 3% maximum brightness. The two plus signs in each panel 
denote the OH maser points. 

SNR 3C 391 was observed with the Chandra ACIS on 03 August 2002. The 
level 1 raw event data were reprocessed to generate a level 2 event file using 
the CIAO software package (version 2.3), resulting a net 60.7 ks exposure. In 
Figures 1 and 2, we show the diffuse X-ray map in 0.3-7 keV and the narrow 
band 1.2-1.5, 1.7-2.0, and 2.3-2.6 keV (including Mg Hea, Si Hea, and S Hea, 
respectively) diffuse emission images. The equivalent width maps (i.e., the 
line-to-continuum ratio maps) for Mg, Si, and S are shown in Figure 3. 
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These X-ray images display the SE-NW elongated morphology, revealing a 
highly clumpy structure of the remnant, with clumps or knots located in both 
the SE and NW parts. Several remarkable, very bright knotty features appear 
on the east and SE border (regions #1 and #2) of 3C 391. These bright knots 
on the border may be small clouds that have recently been shocked by the 
supernova blast wave. A bright enhancement is peaked near the NW border 
(region #3). A complex mixture of knots is seen in the SE part of the remnant 
interior, including at least four bright enhancements indicated as regions #6, 
#7, #8, and #9. 



In Fig.l, arc- or shell-like structures are seen along the northeastern and north- 
ern rim (regions #10 and #11). An X-ray brightened slab at the west rim 
appears to be very close to (just slightly behind) the radio peak emission 
(region #12), and may be related to a small dense region there. On the south- 
west, faint diffuse emission seems to extend out of the radio border. On the 
corresponding location in the narrow-band and EW images of Si and S lines 
(Figs. 2b, 2c, 3b, and 3c), there is a finger-like feature protruding radially out 
of the radio border. This protrusion looks somewhat similar to the apparent 
protrusion in Si and S lines that was recently revealed on the northeastern 
border of Cas A and explained as one of the jets of ejecta (Hwang, Holt, & 
Petre 2000; Hwang et al. 2004). However, the small number of counts collected 
for this feature makes it difficult to determine the metal abundances. 



With the point-like sources 
removed, X-ray spectra were 
extracted from 17 regions 
shown in Fig.l. Most of the 
small-scale regions are chosen 
to include the small features 
of X-ray enhancement such 
as the knots and the faint 
shell like structures. The area 
on the S3 chip outside region 
mtr<| m(h i |u , VXK| ]ii(Kl( , L #17 was used for background. 

The distinct line features, Mg Hea (~ 1.35 keV), Si Hea (~ 1.85 keV), and 
S Hea (~ 2.46 keV) in the spectra (see Fig. 4) indicate the thermal origin of the 
diffuse gas. We find that the spectra of the diffuse gas can be best described 
by the VNEI model with the correction of interstellar absorption. The spectral 
fit results are tabulated in the Table. 




channel energy (keV) 

Fig. 4 Chandra ACIS spectra of the entire remnant of 3C 391 



The absorption column density is found to generally increase across the rem- 
nant from SE to NW, consistent with the existence of a molecular cloud to 
the NW (Wilner et al. 1998). The spectral fits show that the diffuse emission 
from various regions have ionization parameters (n e t) close to or higher than 
10 12 cm™ 3 s. This implies that the hot plasma in the SNR is very close to, or 
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is basically in, the ionization equilibrium. The spectral fits also show that the 
diffuse emission from various regions can be well fitted with solar abundances 
or abundances very close to solar values. The temperature of the gas interior 
to the SNR is generally ~0.5-0.6keV, with only small fluctuations. The gas 
density of each defined region was roughly estimated (as also listed in the 
Table). Apart from the compact knots on the SE and eastern boundary, most 
of the bright knots have a gas density ~ 5-7/ _1 / 2 rf 8 cm" 3 and most of the 
regions along the remnant border have a density ~ l-?>f~ l l 2 d s cm" 3 (where 
d 8 = d/8kpc denotes the distance). The X-ray luminosity in 0.5-10 keV of 
the remnant is ~ 3.5 x lO^c^ergss" 1 . 



VNEI fitting results with the 90% confidence ranges and density estimates 
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a: / is the filling factor of the hot gas. 

b: Making abundances of Mg, Si, and S free parameters apparently improves the fit. 



3 Discussion on the composite appearance 

Here we briefly compare the four mechanisms mentioned in §1 with the prop- 
erties found from our spatially-resolved spectral analysis. 
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i) Projection Effect The molecular cloud is located in the NW but the X- 
ray emission is enhanced not only in the NW half, but also in the SE half. 
Additionally, the variation of the hydrogen column density across the remnant 
implies that the density gradient of the ambient medium seems to be close 
to the projection plane. Therefore, the projection effect does not match the 
observational properties, ii) Radiative Rim The filamentary near-infrared [Fe 
II] and the mid-infrared 12-18 fim [Ne II] and [Ne III] emission along the NW 
radio shell provide some evidence for a radiative cooling at the rim (Reach, 
Rho, & Jarrett 2002). Along the NW border, the shell formation time is t shc ii ~ 
3.8 x lO 3 -^ 14 yr, comparable to remnant's age t = (2/5)(r s /v s ) ~ 4 x 10 3 yr. 
However the X-ray emission along the border (e.g. regions #10, #11, #12, and 
#16), which arises from hot gas (~ 7 x 10 6 K), indicates that a considerable 
amount of gas at the blast shock has not yet suffered significant radiative 
cooling. Additionally, the rim cooling mechanism can not explain what the 
central bright clumpy emission is. iii) Thermal Conduction For 3C 391, the 
thermal conduction scenario is favored by the radiative filaments along the NW 
border and the high ionization timescale implicative of relatively little newly 
shocked material. The conduction timescale t C ond ~ 5.2 x 10 3 (£/2.3 pc) 2 yr 
for the spatial scale £ of order the separation between clumps (typified by 
1' ~ 2.3 pc) is comparable to the remnant's age and hence implies a role 
of conduction in smoothing the interior temperature profile. However, the 
thermal conduction scenario asks for a decrease of temperature and an increase 
of gas density with remnant radius. This predicted behavior is inconsistent 
with the uniform distribution of the gas temperature and density, iv) Cloudlet 
Evaporation The relatively uniform distribution of temperature (even with 
slightly lower values at the center) is actually expected by the White & Long 
(1991) cloudlet evaporation model for model parameters r — >• oo and C/r ^ 3, 
where r is the ratio of the cloud evaporation timescale to the SNR's age 
and C is the ratio of the mass in the cloudlets to the mass of ICM. The 
ratio between the mean density (~ 2f~ 1 / 2 dg 1 ^ 2 cm -3 ) and the density along 
the border (~ 1-3 f^^dg cm -3 ) is basically consistent with that predicted 
in the evaporation model. In fact, the highly clumpy structure unveiled in 
this observation lends support to the conjecture that the ambient molecular 
cloud is inhomogeneous. A combination of very dense clumps and moderately 
dense gas are directly observed in the millimeter molecular lines (Reach & 
Rho 1999). The cloudlets engulfed by the supernova blast wave can act as a 
large reservoir of interior gas by gradual evaporation. This could also explain 
why the interior gas density is much lower than that of the ambient cloud gas 
(~ 30-100 cm -3 ). With this model, the supernova explosion energy is E ~ 0.3- 
1.4 x 10 51 G?!(n /0.3 cm -3 ) ergs. The main difficulty with the model is the high 
ionization age (> 10 12 s) in some regions (e.g. region #1) opposed to the low 
age of the newly evaporated gas. 
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